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Abstract In this review, we describe how high resolu¬ 
tion near infrared spectroscopy and spectro-astrometry 
have been used to study the disks around Herbig Ae/Be 
stars. We show how these tools can be used to iden¬ 
tify signposts of planet formation and elucidate the 
mechanism by which Herbig Ae/Be stars accrete. We 
also highlight some of the artifacts that can complicate 
the interpretation of spectro-astrometric measurements 
and discuss best practices for mitigating these effects. 
We conclude with a brief discussion of the value of long 
term monitoring of these systems. 

Keywords stars: circumstellar matter, individual 
(HD 100546), protoplanetary disks 

1 Introduction 

Herbig Ae/Be stars are intermediate mass (2-8Mq) pre- 
main sequence and zero age-main sequence stars whose 
spectral energy distribution reveals an infrared excess 
(Herbig 1960; Waters & Waelkens 1998; Vieira et al. 
2003). They are the precursors to main sequence A- 
stars many of which have debris disks (Su et al. 2005; 
2006). As the higher mass analogs to classical T Tauri 
stars, Herbig Ae/Be stars provide an important link 
between their lower mass solar analogs and high mass 
stars. 

A great deal of progress has been made in the 
study of the circumstellar environment of Herbig Ae/Be 
stars since the last international conference on Herbig 
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Ae/Be stars (The et al. 1994a). Over the past two 
decades, advances in coronagraphic imaging (Grady et 
al. 1999, 2005, 2014; Weinberger, et al. 1999), the 
(sub)millimeter (Mannings & Sargent 1997; Dent et al. 
2005) and infrared observations (e.g. Meeus et al. 2001, 
2012), and SED modeling (Dullemond et al. 2007) 
have advanced our understanding of this environment. 
While the the geometry of the dust giving rise to the 
infrared excess was controversial two decades ago (The 
et al. 1994b), it is now clear that Herbig Ae/Be stars 
possess flared disks in a Keplerian orbit - usually with 
a puffed up inner rim. Additionally, the study of post- 
main sequence stars that have evolved from A stars in¬ 
dicate that these disks form giant planets perhaps even 
more prodigiously than their lower mass counterparts 
(Johnson et al. 2010). 

In addition to these tools, near infrared spectroscopy 
provides valuable insight to the gas in the inner region 
of the disk. For example, while this region (<50 AU) 
remains difficult to study by direct imaging, the kine¬ 
matic profile of resolved emission lines can serve as a 
surrogate for high resolution imaging by providing spa¬ 
tial information about the gas on milliarcsecond scales 
assuming the gas is in a Keplerian orbit and that the 
stellar mass and disk inclination are known. Ground 
based high resolution (R>20,000) spectroscopy of ro- 
vibrational CO lines have been used to measure the 
distribution of molecular gas relative to dust in the in¬ 
ner disk. This technique was pioneered in the study of 
classical T Tauri stars and embedded YSOs (e.g. Na¬ 
jita et al. 2000, 2007) and has been applied extensively 
to gas in disks around Herbig Ae/Be Stars (e.g. Brit¬ 
tain & Rettig 2002; Brittain et al. 2003, 2007; Blake 
& Boogert 2004; Goto et al. 2006; van der Plas et al. 
2009, 2010; Salyk et al. 2009; Brown, J. et al. 2013). 
The study of ro-vibrational CO lines has also comple¬ 
mented FIR studies of disks around Herbig Ae/Be stars 
by adding parameter space to constrain thermochemi- 
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cal models of disks (e.g. Carmona et al. 2014; Thi et 
al. 2014). 

Spectroscopy of the infrared H I recombination lines 
has also advanced our understanding of Herbig Ae/Be 
stars. Moderate resolution spectroscopy combined with 
interferometry has constrained the size of the emitting 
region of H I to typically within ~0.1 AU (Eisner et al. 
2009; Kraus et al. 2008). A great deal of work has also 
been done to calibrate these lines against the accretion 
luminosity so that they can be used to measure the 
accretion rate of Herbig Ae/Bes (van den Ancker et al. 
2005; Garcia Lopez et al. 2006, Donehew & Brittain 
2011; Mendigutia et al. 2011; Pogodin et al. 2012). 

A particularly exciting development in the use of 
NIR spectroscopy to study gas in disks is the applica¬ 
tion of spectro-astrometry. Spectro-astrometry is the 
measurement of the spatial center of the point spread 
function (PSF) of a spectrum as a function of wave¬ 
length (Bailey et al. 1998a). Because variations in the 
center of the PSF can be measured on much smaller 
angular scales than the diffraction limit of an instru¬ 
ment spectro-astrometry can provide spatial informa¬ 
tion on sub-milliarcsecond scales. This technique was 
initially developed with the advent of CCDs in the 
1980s (see Bailey 1998b for a review) and has since 
been applied extensively to extract additional spatial 
information from spectra of a number of astrophysical 
phenomena such as binaries (Bailey et al. 1998a; Porter 
et al. 2004), YSO jets (Takami et al. 2003, Whelan et 
al. 2004; Davies et al. 2010), and circumstellar disks 
(Acke & van den Ancker 2006; Pontoppidan et al. 2008; 
Brittain et al. 2013; Brown, L. et al. 2013). 

In this review we discuss how modeling line profiles 
spectrally and spatially provide information about the 
distribution of gas in disks around Herbig Ae/Be stars. 
We then highlight two example applications: identifica¬ 
tion of signposts of planet formation and the study of 
the origin of H I emission lines in Herbig Ae/Be stars. 
We follow with a discussion of sources of potential arti¬ 
facts that may complicate the interpretation of spectro- 
astrometry measurements. We show the scope of these 
artifacts and how to acquire the data so as to mitigate 
their effect. We conclude with a brief discussion of fu¬ 
ture prospects for ground based study of gas in disks. 

2 Technique 

2.1 Line Profile Modeling 

Planet forming disks begin their lives with high optical 
depths at most wavelengths. For the minimum mass so¬ 
lar nebula, the surface density at 1 AU is 1700 g cm -2 


or about N(H)=10 27 cm -2 . For a gas to dust ratio and 
grain size distribution representative of the local inter¬ 
stellar medium, only the top 10 -5 of the disk is opti¬ 
cally thin at 5/nm. Grain growth and settling increase 
the depth of the optically thin portion of the disk at¬ 
mosphere, but even so NIR observations only probe a 
thin layer of the disk. The line emission observed in 
the disk arises from this thin warm layer of the atmo¬ 
sphere. Thus observing H 2 can be very challenging be¬ 
cause of the weak oscillator strengths of the quadrupole 
transitions (Bitner et al. 2007; Carmona et al. 2011) 
and large H 2 column densities are needed for detectable 
emission. CO has proven to be the most commonly ob¬ 
served species in disks despite the fact that its abun¬ 
dance is 10 -4 10~ 6 that of H 2 . The ro-vibrational 

transitions of CO are a useful tracer of gas for a num¬ 
ber of reasons. First, the relatively low sublimation 
temperature and high dissociation temperature makes 
it abundant throughout a circumstellar disk. Second, 
even relatively small columns of CO are self shielding. 
Ro-vibrational transitions are also readily observable 
from the ground (moderate Doppler shifts are neces¬ 
sary for the J<11 v=l-0 transitions). 

We can trace the distribution of this warm layer of 
gas in the disk by modeling the kinematic profile of 
the lines. If we assume the gas follows a circular Ke- 
plerian orbit and know the disk inclination and stellar 
mass, then the line profile can be readily calculated (e.g. 
Srnak 1981; Figure 1). If the line is resolved, the inner 
and outer edge of the gas emission can be calculated. 
The half width at zero intensity (HWZI) of the emission 
line corresponds to the projected velocity of the gas at 
the inner edge of the detected region. The half-width 
of the peak separation corresponds to the velocity of 
outer extent of the emitting region. 

Generally the data are modeled by assuming the gas 
distribution is axisymmetric and follows circular Kep- 
lerian orbits, however, this may not accurately describe 
gas in circumstellar disks. For example, the gas in orbit 
may be elliptical giving rise to asymmetric line profiles 
(Regaly et al. 2010) or the gas may arise in a non- 
Keplerian orbit such as a disk wind (Carr 2007; Pon¬ 
toppidan et al. 2011). 

2.2 Spectroastrometry 

Using a relatively new technique, spectroastrometry, 
one can test the assumption that the gas arises in a Ke- 
plerian orbit. Spectro-astrometry has been described in 
a number of papers (e.g. Bailey 1998ab; Porter et al. 
2004; Brannigan et al. 2006; Whelan & Garcia 2008), 
so we limit our discussion to a brief overview of this 
technique - particularly as it applies to circumstellar 
disks. 
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Fig. 1 Schematic of a rotating disk, emergent spectral line 
and spectro- astrometric signal. In the top panel the pro¬ 
jected isovelocity contours are colored. The corresponding 
parts of the emission line are given the same color. The 
half width at zero flux corresponds to the velocity of the 
gas at the inner edge of the disk. The half-width of the 
peak separation corresponds to the velocity of outer extent 
of the emitting region. When the slit of the spectrograph is 
aligned with the semi-major axis of the disk (the rectangle 
superimposed on the schematic), the blue shifted gas is spa¬ 
tially offset toward one end of the disk and the red-shifted 
gas is spatially offset toward the other end of the slit. The 
magnitude of the offset depends on the spatial extent of the 
gas, the radial distribution of the disk flux, and the relative 
brightness of the continuum. 


is not a direct measure of the extent of the gas, rather it 
is the projection of the center of light of a given velocity 
channel along the axis of the slit. 

Typically, the continuum is much brighter than the 
emission lines, so the continuum dominates the center 
of light. The extent of the offset depends on the an¬ 
gular extent of the emission, the asymmetry of the dis¬ 
tribution of the emission along the axis of the slit, the 
brightness of the emission relative to the continuum, 
and the relative position angle of the disk and the slit. 
For example, if the slit is aligned with the semi-major 
axis of an inclined axisymmetric disk, gas in a Keple- 
rian orbit will show emission lines where the red shifted 
side of the line is spatially offset toward one end of the 
slit while the blue shifted side of the emission line is off¬ 
set toward the opposite end of the slit (figure 1). If the 
slit is rotated so that it is aligned with the semi-minor 
axis of the disk then the center of light of the red and 
blue shifted gas will be aligned with the center of the 
PSF and no offset will be measured. 

This technique has provided valuable information 
about gas in disks. For example, high resolution spectra 
of the [O I] A 6300 A lines have been analyzed spectro- 
astrometrically (Acke & van den Ancker 2006). These 
authors show that the line is consistent with gas in Kep- 
lerian orbit about Herbig Ae/Be stars rather than aris¬ 
ing in an outflow as in the case of T Tauri stars. In 
contrast, spectro-astrometric analysis of high resolution 
spectra of ro-vibrational CO emission from strongly ac¬ 
creting YSOs shows that the gas is less extended than 
the narrow profiles indicate if the gas were in a Keple- 
rian orbit (Pontoppidan et al. 2011; Bast et al. 2012). 

In the section that follows we describe two appli¬ 
cations of spectral profile fitting and SA to elucidate 
aspects of disks around Herbig Ae/Be stars. 


While the spatial resolution of a PSF is limited to 
~1.2D/A, the centroid of the PSF can be measured to a 
small fraction of that. For a well sampled Gaussian PSF 
dominated by photon noise, the center of the Gaussian 
can be measured to an accuracy of, 
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0.4 


FWHM 

SNR 


(1) 


where b is the center of the Gaussian, FWHM is the 
full width at half maximum power of the Gaussian, and 
SNR is the signal to noise ratio of the PSF (Brannigan 
et al. 2006). With a SNR of 100 and an AO corrected 
beam of 0".12, the centroid of the PSF can be measured 
to an accuracy of 0.5 milliarcseconds assuming the PSF 
is adequately sampled. The spectro-astrometric signal 


3 Applications 

3.1 Signposts of Planet Formation 

The identification of gas giant planets forming in cir- 
cumstellar disks will represent a major milestone in the 
study of planet formation. By identifying forming plan¬ 
ets we can connect the initial conditions in disks to the 
kinds of planets that form. We can also determine the 
disk radii at which planets form for comparison with 
the orbital radii of planets observed in more mature 
systems. Such observations can also test theoretical 
predictions that planets clear gaps in disks and accrete 
mass through circumplanetary disks, processes that af¬ 
fect the evolution of planetary masses. This approach 
is complementary to the direct detection of forming gas 
giant planets. 
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Direct detection is challenging due to the small an¬ 
gular separation between the forming gas giant planet 
and the host star as well as the high contrast between 
the forming planet and the disk. However, several ex¬ 
citing candidate forming objects have been identified 
by direct imaging (e.g., T Cha Huelamo et al. 2011; 
LkCa 15 Kraus & Ireland 2012; HD 100546 Quanz et 
al. 2013). 



Fig. 2 Synthetic line profile from gas arising from an ec¬ 
centric annulus. The line profile of a line arising from an 
annulus is plotted. The asymmetry of the line changes with 
the eccentricity of the disk. The asymmetry arises from two 
effects. First, the gas near periastron is moving faster than 
gas near apastron. Second, the gas at periastron is closer to 
the star and thus brighter than the more distant gas. Even 
if the line is not resolved, this effect can be inferred from 
the shift of the centroid of the line relative to the systemic 
velocity of the system. 

Another approach to identifying forming gas giant 
planets is to identify dynamical markers of tidal inter¬ 
actions between the disk and forming object. Hydro- 
dynamic modeling of planet disk interaction indicates 
that a planet with a mass greater than ~3Mj up i ter on 
a circular orbit can induce an eccentricity in the disk 
in the vicinity of its orbit (e.g. Kley & Dirksen 2006). 
Simulations show that the eccentricity can be as high 
as 0.25 and falls off as approximately r~ 2 . The profile 
of a line arising from an eccentric annulus will gen¬ 
erally not be symmetric. The gas at periastron will 
be faster than the gas at apastron and generally the 
temperatures ■will differ as well. This results in a dis¬ 
tinctive asymmetric line profile (figure 2; see Regaly et 
al. 2010 and Liskowsky et al. 2012). Interestingly, the 
semi-major axis of the disk precesses very slowly (~T0° 
per 1000 orbits), so the line profile should remain con¬ 
stant over decadal timescales. Such a line profile has 
been observed in OH spectra from the Herbig Ae bi¬ 
nary V380 Ori (Fedele et al. 2011) and the Herbig Be 


star HD 100546 (Liskowsky et al. 2012). In the case of 
V380 Ori, it is likely a stellar companion that is driving 
the eccentricity whereas in the case of HD 100546 it is 
likely that a substellar object is driving the eccentricity. 
In principle, an asymmetric line profile could also arise 
from the uneven sampling of a resolved axisymmetric 
disk (e.g. Hein Bertelsen et al. 2014) which we will 
discuss in section 3. 

Another way to identify ongoing planet formation 
is by detecting emission arising from a circumplane- 
tary disk feeding the planet. Circumplanetary disks 
can be discovered by their line emission (Brittain et 
al. 2013; 2014) and also potentially by their dust con¬ 
tinuum emission (Zhu 2014; Isella et al. 2014). Sim¬ 
ulations of gas giant formation indicate that the disk 
feeding the gas giant planet will have a radius roughly 
one-third of the Hill Sphere (Quillen & Trilling 1998; 
Ayliffe & Bate 2009ab; Martin & Lubow 2011), though 
the full accretion envelope is likely intrinsically three 
dimensional with complicated flow patterns (e.g. Tani- 
gawa et al. 2012; Ayliffe & Bate 2012; Gressel et al. 
2013). Thus a 5Mj up it er planet forming at an orbital 
radius of 10 AU should have a circumplanetary disk 
with a radius of 0.3 AU - comparable to the size of the 
emitting region of ro-vibrational CO emission observed 
toward classical T Tauri stars (Najita et al. 2003). 

In the case of Herbig Ae/Be stars with transition 
disks, hot bands of ro-vibrational CO emission are ob¬ 
served from R>10 AU (e.g. HD 179213, van der Plas 
2010; HD 100546, Brittain et al. 2009; HD 141569 Brit¬ 
tain et al. 2003; HD 97048, van der Plas et al. 2009; 
and Oph IRS 48, Brown et al. 2012) due to the UV 
fluorescence of the gas. To detect the IR ro-vibrational 
emission lines that result from this process, a large area 
is needed because these transitions are exceeding opti¬ 
cally thin. For example, the line flux of CO emission 
arising from HD 100546 is comparable to that observed 
from classical T Tauri stars in Taurus (compare for ex¬ 
ample, Brittain et al. 2009 and Najita et al. 2003). 
However, the emitting region of the CO observed to¬ 
ward HD 100546 extends over 50 AU while the emis¬ 
sion from the classical T Tauri stars typically originates 
within 1 AU of the star. Thus UV fluoresced gas orig¬ 
inating from a disk with a radius less than 1 AU will 
reveal ro-vibrational CO lines with luminosities 3-4 or¬ 
ders of magnitude fainter than what is observed from 
the transition disks mentioned above. 

The contrast between the gas conditions which pro¬ 
duce UV fluorescence (broadly distributed, perhaps 
tenuous, gas) and collisionally excited emission (warm, 
dense material) provides a way to identify circumplan¬ 
etary disks. The surface density and temperature of 
the circumplanetary disk is expected to be enhanced 
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Fig. 3 Evolution of CO emission line profile and spectro-astrometric signal (Brittain et al. 2014). In the first column 
the v=l-0 P26 CO line observed in 2003 is plotted (black) and compared to the same line observed in 2006 (red), 2010 
(green), and 2013 (cyan). The Doppler shift of the excess emission can be seen in this column. In the second column, 
the spectro-astrometric signal is plotted for each epoch along with the model (dot-dashed line) and the residual. The 
spectro-astrometric signal is modeled by adding an extra source of emission in the inner edge of the disk at the position 
corresponding to the Doppler shift of the excess (illustrated in the third column). While the inner edge of the disk is roughly 
130mas from the star, the center of light of the CO emission is offset from the center of the PSF ±15mas along the axis of 
the slit. The extra source of emission in orbit around the star causes the center of light to shift. This source of emission 
was hidden by the circumstellar disk in 2003. 


relative to the feeding material (Klahr & Kley 2006; 
Matin & Lubow 2011). If the gas is of order 10 3 K and 
optically thick, then the flux of the emission from the 
circumplanetary disk will be similar to what is observed 
from classical T Tauri stars. 


A Herbig Ae/Be star with a transition disk created 
by a forming gas giant planet should then give rise to 
v=l-0 lines that are a blend of emission from the cir¬ 
cumstellar disk and circumplanetary disk while the hot- 
band transitions will be overwhelmingly dominated by 
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the UV fluoresced gas in the circumstcllar disk. Thus 
as the circumplanetary disk orbits the star, the shape 
of the v=l-0 lines will change with the motion of the 
circumplanetary disk. The hotband lines provide a con¬ 
stant fiducial profile against which to measure the evo¬ 
lution of the v=l-0 lines. This behavior, variable v=l 
0 profiles with constant hotband profiles, is observed in 
HD 100546 (Brittain et al. 2013; 2014). 

In January 2006 the red side of the v=l-0 P26 line at 
0-10 km s -1 brightened relative to what was observed 
in January 2003 (second row of figure 3). The Doppler 
shift of this excess emission was 6±1 km s _1 . In De¬ 
cember 2010, the v=l-0 P26 was still brighter than it 
was in January 2003, however, the Doppler shift was 
1±1 km s -1 (third row of figure 3). In March 2013 
the excess emission was Doppler shifted -6±1 km s -1 . 
Assuming the source of the excess emission is in a Kep- 
lerian orbit, the projected Doppler shifts and epochs of 
the observations place the emission near the inner edge 
of the outer disk (R~13 AU). One way to test whether 
the shift of the excess v=l-0 line emission is due to gas 
in a Keplerian orbit and determine the radial position 
of the emission is to use spectro-astrometry. 

Consider an inclined axisymmetric disk of gas in a 
circular Keplerian orbit. If the slit is aligned near the 
semi-major axis of the inclined disk, then the blue side 
of the emission lines formed in the disk will be offset in 
one direction along the slit axis and vice versa for the 
red side of the line, and the offsets would be symmetric 
(figure 1). This was observed for the v=l-0 P26 line in 
2003 (top row of figure 3). But what happens if we add 
an extra localized source of emission at a given velocity 
to our disk? In this case, the center of light at the 
velocity of the extra source will be shifted closer to the 
added source of emission. 

The extent and direction of the shift will depend on 
the location of the center of light along the slit axis. 
In the case of HD 100546, the first three observations 
were taken with the slit in its default position angle of 
90° (Brittain et al. 2013). The position angle of the 
disk is 140° (Ardila et al. 2007). In 2006 not only was 
the red side of the line brighter relative to its bright¬ 
ness in 2003, the center of the PSF of the red side of 
the line was less extended in 2006 than it was in 2003 
(second row of figure 3). As it turns out, for the com¬ 
bination of the slit PA and disk PA, an excess source 
of emission that orbits close to the inner disk edge of 
the outer disk (R^13 AU) with a projected velocity of 
6 km s -1 will be located near the center of the slit, 
just as required to explain the spectro-astrometric sig¬ 
nal (second row of figure 3). In 2010, the spectral pro¬ 
file of the extended emission was broader and shifted to 
—1±1 km s -1 . Placing this emission in the disk near the 


inner rim at the orbital position corresponding to this 
shift reproduces the spectro-astrometric signal (third 
row of figure 3). The same is observed for the observa¬ 
tions acquired in 2013. From 2006 to 2010 to 2013, the 
Doppler shift of the emission, location of the emission 
projected along the slit axis, and period of the emission 
are all consistent with a source of thermal CO emission 
in a circular Keplerian orbit near 13 AU (Brittain et al. 
2014). 

3.2 The Origin of H I lines 

A second important application of spectro-astrometry 
to the study of Herbig Ae/Be stars is the study of the 
origin of the H I recombination lines. The origin of 
these lines remains uncertain (Kraus et al. 2008): do 
they form in the stellar wind, disk wind, decretion disk, 
accretion flow, or some combination thereof? If the 
H I lines are connected to accretion, then they may 
serve as a useful proxy of the stellar accretion rate (e.g. 
Garcia Lopez et al. 2006). Several studies have at¬ 
tempted to calibrate various emission line diagnostics 
to the accretion luminosity (e.g. Donehew & Brittain 
2011; Mendigutfa et al. 2011, 2012, 2013; Oudmaijer 
et al. 2011; Pogodin et al. 2012; see also Salyk et al. 
2013). Typically, these studies convert the veiling of 
the Balmer discontinuity to an accretion luminosity us¬ 
ing a magnetospheric accretion model (e.g. Muzerolle 
et al. 2004) and compare this to the luminosity of var¬ 
ious emission lines. For example, Donehew & Brittain 
(2011) apply this method and find that the relation¬ 
ship between the luminosity of Bry and the accretion 
luminosity determined for Herbig Ae stars is consistent 
with the relationship found for classical T Tauri stars 
(Muzerolle et al. 1998) and intermediate mass T Tauri 
stars (Calvet et al. 2004); however, the Herbig Be stars 
do not follow the same trend. 

Whether the correlation between the luminosity of 
the H I emission lines and the inferred accretion lumi¬ 
nosity is physical or fortuitous remains to be seen. One 
way to investigate this issue is to measure the vari¬ 
ability of the veiling of the Balmer discontinuity and 
emission line luminosities simultaneously (Mendigutia 
et al. 2013). With the advent of broad coverage 
spectrographs such as X-Shooter on the VLT (cover¬ 
ing ~ 0.3 — 2.4/zm), it is now possible to simultane¬ 
ously observe the veiling of the Balmer discontinuity 
and the emission line diagnostics that fall within this 
range (Oudmaijer et al. 2011; Mendigutia et al. 2014). 
Ongoing studies along these lines will be able to deter¬ 
mine whether or not the veiling and line emission are 
connected. However, they will not be able to establish 
that Herbig Ae/Be stars accrete magnetospherically - 
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Fig. 4 Spectro-astrometric measurement of Pa (3 (left column) and Br 7 (right column) with NIFS on Gemini North. 
The spectral profile of the line is plotted (black) above the spectro-astrometric signals. The red lines correspond to the 
centroid measurements made along the slitlet (the x-direction) and the blue lines refer to the centroid measurements made 
perpendicular to the slitlet (the y-direction). No significant signal is detected for the HAe stars, but a signal is detected for 
the HBe stars. Modeling the line profiles and signals is an ongoing project (Adams et al. 2015 - in prep.). 


as assumed in the calculation of the accretion luminos¬ 
ity - or identify where the H I emission lines originate. 

At first blush the assumption that Herbig Ae stars 
accrete magnetospherically would seem unjustified. 
Herbig Ae/Be stars are not fully convective thus one 
would not expect that they would generate the strong 
kilogauss magnetic fields characteristic of their lower 
mass analogs (Johns-Krull 2007). Indeed, studies of 
the magnetic properties of Herbig Ae stars indicate 


that they do not possess strong well ordered surface 
fields (Alecian et al. 2014). On the other hand red 
shifted hydrogen lines are observed in Herbig Ae stars 
with velocities of several hundred km s^ 1 indicative 
of material in free fall onto the star (the tell tale sig¬ 
nature of magnetospheric accretion; Guimaraes et al. 
2006). Thus the assumption that at least some Herbig 
Ae/Be stars accrete magnetospherically is not without 
warrant. 
















































Before the recombination H I lines can be used re¬ 
liably as a proxy for the stellar accretion rate, it is 
crucial that we understand the physical origin of these 
lines. One promising development is the use of spectro- 
astrometric measurements of H I emission using in¬ 
tegral field spectrographs. This method has recently 
been applied to massive YSOs. Davies et al. (2010) 
measured Bry emission from the massive YSO W33A 
with NIFS on Gemini North. They were able to mea¬ 
sure the centroid of the PSF with an accuracy of 0.1 
milli&rcseconds and show that the line originates in an 
outflow. 

Application of the same method to Herbig Ae/Be 
stars is beginning to show promising results. For ex¬ 
ample (Adams et al. 2015), show they can achieve a 


root mean square fidelity of 100/iarcseconds with their 
observation of Herbig Ae/Be stars. This is within a fac¬ 
tor of a few of being able to distinguish emission from 
a magnetosphere, outflow, and gas in a Keplerian orbit 
in the inner 0.1AU. A study of four Herbig Ae/Be stars 
finds SA signatures from H I from the HBe stars but 
not the HAe stars in the sample (Figure 4). By combin¬ 
ing these data we are also able to show that our fidelity 
is limited by photon noise. Increased integration time 
should provide a means to improve the sensitivity of 
these measurements by a factor of a few and thus dis¬ 
criminate between a wind, disk, or funnel flow origin 
for the nearest Herbig Ae/Be stars. Making such high 
precision measurements requires attention to potential 
artifacts in the data. 
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Fig. 5 Illustration of artifacts resulting from pointing errors. In row A, the isovelocity contours are plotted in red. Also 
plotted are the same contours convolved with two dimensional Gaussian profile with a full width at half max equal to 0A6. 
The portion of the disk sampled by the slit on the Phoenix spectrograph is marked with green dot-dashed lines. In row B, 
the resultant line profile arising from the disk is plotted. In the first panel, the PSF of the disk is centered in the slit. In the 
next the center of the PSF is offset by one pixel, and so on. In row C, the spectro-astrometric measurement of the PSF is 
plotted. Each panel corresponds to the offset noted in the row above. When the disk emission is resolved and the PSF is not 
properly centered in the slit, the profile is asymmetric. This has a minor effect on the line profile and spectro-astrometric 
signal unless the offset is substantial. Misaligning the center of the PSF and slit results in about a 20% error when offset is 
half a slit width and 50% error when the PSF is offset a full slit width. With care the center of the PSF can be measured 
and placed in the center of the slit with an accuracy of about 0.5 pix. 
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4 Pitfalls 

Several authors have summarized best practices for us¬ 
ing SA to study phenomena such as binaries and out¬ 
flows (Bailey et al. 1998; Baines et al. 2004; Porter 
et al. 2004; Brannigan et al. 2006; Whelan & Garcia 
2008). Because SA involves making high precision mea¬ 
surements of the center of the PSF, even small artifacts 
can lead to spurious signals. We summarize the prin¬ 
cipal best practices extracted from the work of these 
authors: 

1. Observe the source twice at each desired PA - once 
at the original position and then with the slit ro¬ 
tated 180°. Instrumental artifacts should rotate 
with the instrument. If they are reproducible, then 
they should cancel when the signals are combined. 

2. Observe a standard star as a secondary check for 
instrumental artifacts 

3. Acquire very high signal to noise ratio flats every 
time the grating is moved. If the grating in the spec¬ 
trograph moves when the telescope slews, it may 
be necessary to retake flats with each observation. 
Small errors in the flat-fielding can shift the center 
of the PSF and result in spurious signals. 

4. The slit width should be narrow compared to the 
PSF. If the PSF does not fill the slit and the PSF is 
not perfectly symmetric, it can result in a spurious 
signal that will not be removed by combining with 
the SA signal acquired with the instrument rotated 
180°. 

5. Acquire spectra over a wide spectral range. It is 
crucial that the continuum be well sampled to ensure 
a baseline measurement is available for comparison 
to the region covered by the spectral line. This can 
prove to be challenging when observing strong H I 
lines from Herbig Ae/Be stars with high resolution 
spectrographs as the the full width at zero intensity 
of the line can be comparable to the spectral grasp 
of the order. 

In what follows, we highlight two additional is¬ 
sues that may complicate the interpretation of spectro- 
astrometric measurements. First we discuss the chal¬ 
lenges of analyzing spectro-astrometric signals from re¬ 
solved sources and second we explore the challenges of 
making SA measurements from spectra acquired with 
an IFS. 

When observing molecular emission from nearby 
disks, it is possible that the emission will be spatially 
resolved. For example, the molecular NIR emission ob¬ 
served from HD 100546 extends from ~10±2 AU to 
beyond 50 AU (van der Plas et al. 2009; Brittain et al. 
2009; Carmona et al. 2011; Liskowsky et al. 2012). In 
addition there is a compact ring of solids that extends 





Fig. 6 Illustration of the effect of an asymmetric PSF. 
In panel A, the PSF of a star acquired with PHOENIX on 
Gemini South is plotted (red contours). A synthetic, elon¬ 
gated profile is also plotted (black contours). An emission 
line was synthesized assuming it arose from gas in a disk 
in a Keplerian circular orbit. In panel B the resultant line 
profiles that resulted from the PSF plotted in panel A are 
plotted. In panel C the spectro-astrometric signal is simi¬ 
larly plotted. While the focus of the PSF has only a minor 
effect on the spectral line profile, the effect on the measure¬ 
ment of the spectro-astrometric signal is significant. Proper 
focusing of the source is important for acquiring reliable 
signals. 
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from 0.25 AU to less than 0.7 AU that dominates the 
continuum in the NIR (Benisty et al. 2010; Mulders et 
al. 2011, 2013; Panic et al. 2014). At the distance to 
HD 100546 (97±4 AU; van Leeuwan 2007), the inner 
hole of the disk subtends 0".28 ± 0".03 along the semi¬ 
major axis. Adopting an inclination of 42° (Ardila et 
al. 2007; Pineda et al. 2014) indicates that the inner 
hole subtends 0".19 ± 0".03 along the minor axis of the 
disk. 

Spectroscopic observations of molecular lines from 
the disk using adaptive optics and 0".2 slit will result 
in the occultation of a significant fraction of the disk. 
Small changes in the positioning of the slit can give rise 
to varying line shapes (Hein Bertelsen et al. 2014). Re¬ 
construction of the line profile representing the under¬ 
lying structure of the disk thus requires precise knowl¬ 
edge of the slit position. Even pointing uncertainties 
of order OAOS can undermine the interpretation of the 
spectrum. In addition, since the continuum is unre¬ 
solved, the PSF of the continuum and the PSF of the 
emission lines is not identical. This can give rise to 
artifacts as well since the PSF does not fill the slit. 

In such cases it can be advantageous to work at lower 
spatial resolution and wider slit widths. Here we com¬ 
pare the effect of slight pointing offsets on the shape of 
the spectral lines to the effects described by Hein Ber¬ 
telsen et al. (2014). We consider conditions with seeing 
of 0".6 and a four pixel slit with a width of 0".34 (the 


slit width of Phoenix on Gemini South; iHinkle et al, 
20031 2000l 1998h . If the slit is offset from the center of 


the PSF, an asymmetric PSF and spectro-astrometric 
signal does result (figure 5); however, to achieve a devi¬ 
ation of 20% from the original measurement, the PSF of 
the continuum must be shifted as much as 2 pixels from 
the center of the slit. With care the star can be cen¬ 
tered on the slit to within ±0.5pixels (about ±0".05 for 
Phoenix on Gemini South and CRIRES on the VLT). 
While such an offset can be significant when using adap¬ 
tive optics (a PSF FWHM- 0".15) and 0".2 slit (Hein 
Bertelsen et al. 2014) it is much less significant when 
the seeing is 0".6. 

An additional concern arises from making measure¬ 
ments of spectro-astrometric signals from an uncor¬ 
rected PSF. For example, an asymmetric PSF can give 
rise to a spurious spectro-astrometric signal. We mod¬ 
eled the effect of a poorly focused beam convolving our 
disk model with an observed PSF and an asymmetric 
PSF (figure 6a). The line profile shows only a modest 
effect due to the occultation of part of the disk (figure 
6b). The spectro-astrometric signal, however, shows a 
significant asymmetry (figure 6c). Perhaps of greater 
concern is that rotating the instrument will not remove 
artifacts arising from this effect if it is not steady (say 


due to variable seeing). Imaging the PSF throughout 
the observations is a helpful way to identify problems 
with a variable PSF. 


5 Conclusions 

Spectroscopy has been used to elucidate many aspects 
of Herbig Ae/Be stars since the previous conference 
on Herbig Ae/Be stars held in 1994; however, there is 
much exciting work left to be done. Over the past two 
decades, we have moved from detecting and character¬ 
izing line emission arising from circumstellar material 
to using these features to elucidate the star-disk connec¬ 
tion and identify signposts of forming gas giant planets. 
The application of spectro-astrometric measurement to 
the study of disks is still in its infancy. While this 
technique shows remarkable promise for revealing sub- 
rm/Harcsecond spatial information from spectra, atten¬ 
tion is needed to avoid several pitfalls that can com¬ 
plicate the interpretation of such data. The detection 
of companions with orbits of order ~10 AU by mea¬ 
suring the Doppler shift of emission arising from the 
circumplanetary disk requires observations on decadal 
timescales. Archival data is thus a crucial resource fa¬ 
cilitating the long term monitoring of these systems. 
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